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Abstract. We study the response of the solar atmosphere 
to excitations by large amplitude acoustic waves with 
radiation damping now included. Monochromatic adia- 
batic waves, due to unbalanced heating, generate continu- 
ously rising chromospheric temperature plateaus in which 
the low frequency resonances quickly die out. All non- 
adiabatic calculations lead to stable mean chromospheric 
temperature distributions determined by shock dissipa- 
tion and radiative cooling. For non-adiabatic monochro- 
matic wave excitation, a critical frequency v cr ~ 1/25 Hz 
is confirmed, which separates domains of different reso- 
nance behaviour. For waves of v < v cr , the resonances 
decay, while for waves of v > v cr persistent resonance 
oscillations occur, which are perpetuated by shock merg- 
ing. Excitation with acoustic frequency spectra produces 
distinct dynamical mean chromosphere models where the 
detailed temperature distributions depend on the shape of 
the assumed spectra. The stochasticity of the spectra and 
the ongoing shock merging lead to a persistent resonance 
behaviour of the atmosphere. The acoustic spectra show 
a distinct shape evolution with height such that at great 
height a pure 3 min band becomes increasingly dominant. 
With our Eulerian code we did not find appreciable mass 
flows at the top boundary. 

Key words: hydrodynamics - shock waves - waves - 
sun: chromosphere - sun: oscillations 



1. Introduction 

A pronounced signal of velocity and temperature fluctua- 
tions in the solar chromosphere, seen in the Ca II H and 
K, Ha and the Ca II infrared triplet lines, are the 3 min 
(y = 5.5 mHz) oscillations. For detailed reviews of the 
3 min oscillations see Deubner (1991), Fleck & Schmitz 



(1991), Rutten & Uitenbroek (1991), Rossi et al. (1992), 
Carlsson & Stein (1994) as well as Rutten (1995, 1996). 

Fleck & Schmitz (1991) were the first to show that 
the 3 min oscillations can be explained as the basic cut- 
off frequency resonance of the chromosphere. This view 
has been confirmed both analytically and numerically by 
Kalkofen et al. (1994) as well as by Sutmann & Ulmschnei- 
der (1995a, b; henceforth called Papers I, II, respectively). 
For extended recent analytical work see also Sutmann et 
al. (1997; Paper III). The atmospheric resonance is due to 
a local oscillation of gas elements around their rest posi- 
tions in hydrostatic equilibrium. 

Previously, Leibacher & Stein (1981) in the wake of the 
very successful interpretation of the photospheric 5 min 
oscillations had explained the 3 min oscillations as a cav- 
ity mode. It is now universally accepted that the 5 min os- 
cillations are acoustic waves trapped in a subphotospheric 
cavity between the temperature minimum and a refracting 
temperature rise in the solar interior. Leibacher & Stein 
suspected that the chromospheric 3 min oscillation might 
be likewise explained by a chromospheric cavity, acting 
between the temperature minimum and the temperature 
rise of the chromosphere-corona transition layer. 

In our previous work, we have argued against a chro- 
mospheric cavity explanation, because the effect and via- 
bility of the 3 min oscillation as a local resonance both 
by the above mentioned analytical work and by time- 
dependent numerical calculations is now well established. 
Nonetheless, it is highly likely that a chromospheric cav- 
ity resonance might occur in addition to the local reso- 
nance. Observationally, the 180° phase jump in the Na I 
D line points to such a cavity (Fleck & Deubner 1989). 
Theoretical work by Fleck & Schmitz (1991), Kalkofen et 
al. (1994), Carlsson & Stein (1994), Papers I & II as well 
as by Cheng & Yi (1996), however, has only been one- 
dimensional (1-D) and often did not even include the ex- 
istence of the transition layer in the analysis. It can easily 
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be seen in 1-D simulations that acoustic waves will not be 
reflected by a steep temperature gradient. But this is only 
the case because of the very special geometry of purely 
vertical wave propagation. As soon as the wave propaga- 
tion is inclined to the vertical, the temperature gradient 
acts to refract the acoustic wave field and the transition 
layer becomes a reflecting boundary. 

A third way to simulate the 3 min oscillation phe- 
nomenon is to avoid dwelling on mechanisms and simply 
excite the solar atmosphere with velocity fluctuations de- 
rived from Doppler shift observations from a low lying Fe 
I line (Lites et al. 1993). Impressive work using this ap- 
proach aimed to explain the chromosphcric bright point 
phenomenon has been carried out by Carlsson & Stein 
(1994). Similar work using the Lites et al. observations but 
employing a different numerical code has recently been re- 
ported by Cheng & Yi (1996). These latter authors argue 
that the power in the 3 min band is already included in 
the observed signal and that high frequency acoustic wave 
power will not appreciably contribute to the 3 min band. 
This claim has to be taken with some caution as the code 
of these authors, different from that of Carlsson & Stein, 
does not have an adaptive mesh capability and thus does 
not allow to treat shocks accurately. Consequently, the 
code is not able to describe basic acoustic wave proper- 
ties like the limiting shock strength behaviour and shock 
merging. Therefore, the code cannot correctly describe the 
transfer process of high frequency wave power into lower 
frequencies via shock formation and shock merging. 

It is precisely this process of shock merging and the 
resulting transfer of high frequency acoustic wave power 
into the low frequency range, which we want to study 
in greater detail in the present work. We continue here 
the work of Paper II. The reason for this interest is that 
in recent years it became clear that the acoustic energy 
generation in late-type stars is strongly tied to the Kol- 
mogorov type turbulent cascade generated by the surface 
convection zones of these stars, which has an extensive 
high frequency tail (Musielak et al. 1994). It is thus of 
fundamental importance to understand the evolution of 
the acoustic wave spectrum with atmospheric height. It 
should be noted, however, that it is inherently difficult 
to observe the high frequency acoustic wave power be- 
cause the line contribution functions extend over several 
pressure scale heights (Judge 1990) and one consequently 
has small modulation transfer functions (e.g. Deubner et 
al. 1988, Ulmschneider 1990). Another important issue is 
the role of the high and low acoustic wave frequencies in 
the chromosphcric heating. Finally there is the problem 
of how the observed prominent 3 min spectral component 
is to be related to a relatively flat acoustic spectrum pro- 
duced by the Kolmogorov type turbulence in the convec- 
tion zone obtained by Musielak et al. (1994). 

It has been recognized for some time (Rammacher & 
Ulmschneider 1992, Fleck & Schmitz 1993, Kalkofen et 
al. 1994) that for adiabatic waves the process of shock 



merging can greatly modify the atmospheric response and 
will eventually lead to a complete transformation of high 
frequency power into low frequency power. In Paper II, 
confirming work by Kalkofen et al. (1994), it was shown 
that the spectrum of adiabatic acoustic waves continu- 
ously changes towards longer wavelengths with increasing 
atmospheric height, and that this is due to the process 
of shock merging by which short-period shock waves are 
converted to long-period waves. In Paper II we found 
that for monochromatic wave excitation, a critical fre- 
quency v cr exists, below which the monochromatic waves 
stay dominant and above which monochromatic waves are 
completely obliterated by the generation of resonances 
through shock merging. Cuntz (1987), calculating time- 
dependent wave models for Arcturus, also considered ra- 
diative damping. He found that short-period waves are 
converted to long-period ones via shock merging. In this 
work, however, fully developed shocks are already inserted 
at the inner boundary of the atmospheric shell and the 
development of photospheric acoustic waves into shocks is 
not calculated explicitly. 

In the present work we ask whether the theoretical 
results of Paper I and II persist when more realistic non- 
adiabatic wave calculations are performed and the time 
span of the calculations are extended. We also want to in- 
crease the spectral resolution to make more accurate pre- 
dictions about the height evolution of the frequency spec- 
tra. Additionally, we want to investigate the behaviour 
of monochromatic waves of much higher frequency In 
Sect. 2 we briefly outline our numerical methods. Section 3 
presents the results and Sect. 4 gives our conclusions. 



2. Wave calculation method, Fourier analysis, at- 
mosphere model 

The method of wave computation has been described in 
Papers I and II and we refer to these papers for further 
details. The time-dependent hydrodynamic equations are 
solved using the method of characteristics. However, dif- 
ferent to Papers I and II we now use an Eulerian descrip- 
tion (Cuntz & Ulmschneider 1988). We follow the develop- 
ment of the originally linear wave (introduced at the bot- 
tom of the atmosphere by a piston) to the point of shock 
formation and beyond. The shocks are treated as discon- 
tinuities and are allowed to grow to arbitrary strength. 
They are also permitted to merge with each other. We 
consider an atmospheric slab with a height of 2220 km 
and use a total of 593 grid points with a grid spacing of 
3.75 km. In addition to the fixed number of regular grid 
points, there is an arbitrary number (typically 30) of shock 
points, which arc allowed to move between the regular grid 
points according to the speed of the shocks. Since we solve 
the hydrodynamic equation in Eulerian form, our compu- 
tational domain does not move with height and artificial 
net outflows are avoided. 
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Fig. 1. Power spectrum (left panel, drawn) and velocity (right panel) at height z — 2000 km for an excitation with a monochro- 
matic adiabatic wave of period P = 15 s and flux Fa = 1 • 10 8 erg cm~ 2 s _1 . The power spectrum at height z = km, scaled 
by 0.1, is shown dotted. The time span for the Fourier analysis is from 500 to 2548 s, sampled by 1 s intervals. 




Fig. 2. Temperatures as function of height at times t = 500, 2000, 3500 s for a monochromatic adiabatic wave with an acoustic 
energy flux of Fa = 1 • 10 8 erg cm~ 2 s _1 and a period of P = 15 s (left panel) and P = 30 s (right panel). 



Different from Papers I and II, we now consider ra- 
diation damping (see also Rammachcr & Ulmschneider 
1992, and Ulmschneider et al. 1992). The radiative trans- 
fer and statistical equilibrium equations are solved for the 
H continua and the Mg II k and Ly a lines. The lines 
are treated considering complete redistribution and using 
an operator splitting method (Buchholz et al. 1994). The 
radiative emission by the Mg II k line is scaled such as 
to simulate a realistic total chromosphcric radiation loss. 
The usage of Fourier analysis and the evaluation of ve- 
locity power spectra is the same as described in Papers I 
and II, and we refer the reader to these papers for details. 
The atmosphere models are discussed in Paper II. For our 
calculations we take a combined H - , Mg II k and Lya 
radiative equilibrium model as initial atmosphere. 



3. Results 

To investigate the properties of monochromatic waves we 
excited our H _ , Mg II k and Lya radiative equilibrium 
model by a piston at the bottom with waves of various 
periods and an acoustic flux of Fa = 1 ■ 10 8 erg cm~ 2 s -1 . 
Compared to Paper II we now calculate our models with 
a 2.5 times higher spatial resolution. As our time step 
is largely determined by the Courant condition, this in- 
creases the temporal resolution by the same factor. 

3.1. Monochromatic adiabatic wave excitation 

Because the wave calculation with P — 15 s in Fig. 2 of 
Paper II was carried out with a low spatial resolution, we 
first calculate the adiabatic case again, using now about 
25 points per wavelength. Figure 1 shows the power spec- 
tra both at the bottom z = km (dotted) and at Euler 
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Fig. 3. Shock merging heights as function of time for 
monochromatic adiabatic waves with an acoustic energy flux of 
Fa = 1-10 8 erg cm~ 2 s _1 and a period of P — 15 s (diamonds) 
as well as P = 30 s (plus signs). 



height z = 2000 km (left panel) together with the velocity 
at the latter height (right panel). The Fourier spectrum 
was computed for the time span 500 to 2548 s with the 
same sampling rate of 1 s as in the models of Paper II. 
Comparison with Fig. 2a of Paper II shows that, while 
that power spectrum is strongly peaked at about 8 mHz, 
our present spectrum has extended power between 3 and 
30 mHz. This indicates that by decreasing the resolution 
a much narrower resonance oscillation is obtained, which 
can be explained by the fact that with higher spatial res- 
olution the shock merging process is more individualistic, 
leading to a more pronounced low frequency spectrum. 

3. 2. Effect of the generated chromosphere model, adiabatic 
case 

In Paper II we did not discuss the important effect that, 
different from the linear waves in Paper I, non-linear 
waves will lead to significant heating, which consider- 
ably modifies the initial atmosphere and produces a dy- 
namical chromosphere model. Figure 2 shows tempera- 
tures as function of height for adiabatic wave calculations 
with periods of P = 15 and 30 s and an energy flux of 
Fa = 1 • 10 8 erg cm~ 2 s _1 . The waves are displayed at 
times t — 500, 2000 and 3500 s. It is seen that shock for- 
mation occurs in both models at heights near z = 400 km. 
Within a short distance above this height, the waves at- 
tain a sawtooth profile with limiting shock strength and 
produce significant heating. 

An interesting feature visible in Fig. 2 is that a chro- 
mospheric temperature plateau is formed which continu- 
ously rises as function of time. In the 15 s wave case the 
mean chromospheric temperature increases at a rate of 
800 K every 500 s, from 5000 K at t = 500 s to 10000 K 
at t = 3500 s. For the 30 s wave the increase is even faster 



with 2500 K per 500 s, from 6000 K at 500 s to 21000 K 
at 3500 s. The reason for this continuously rising chromo- 
spheric temperature plateau is that despite of a perpet- 
ual heat input, the atmosphere has no way of cooling as 
the waves are treated adiabatically. At a given height in 
the chromosphere each passing shock deposits an entropy 
jump AS per wave period which continuously increases 
the entropy. As for a given mass element the gas pressure 
stays roughly constant, the entropy increase leads directly 
to a temperature increase in the element. 

Let us now discuss the effect of the atmosphere on 
the resonance oscillations. As discussed in Paper II and 
as seen in Fig. 1 the initial switch-on effect dies out after 
about 1200 s at the height of 2000 km. Yet in Fig. 1 and 
in the non- adiabatic cases in Fig. 4 with periods P = 15, 
20, 30 s, there is a lot of low frequency power at times 
t > 1200 s. This is due to the generation of resonances 
by shock merging events. Note that, as discussed in Paper 
II and by Rammacher & Ulmschneider (1992), resonance 
oscillations are both the cause and the result of shock 
merging events. This is because shocks, riding on the high 
temperature part of the low frequency resonant wave pro- 
file, catch up with shocks propagating with lower speed on 
the low temperature part of that wave. The shock merging 
events kick on new transients, which leads to a prolonga- 
tion of the resonances to times far beyond the decay time 
of the switch-on effect. 

We find that the atmospheric model strongly affects 
the shock merging properties. The continuously increasing 
chromospheric temperature plateau causes shock merging 
to occur at progressively greater heights. This is because 
(see Fig. 2) when two shocks of relatively similar strength 
move up a steep temperature rise, the second shock can- 
not easily catch up with the first as it moves with greater 
speed. Only when the temperature plateau is reached, can 
the second shock catch up. The persistently growing tem- 
perature of the plateau leads to an ever increasing height 
where the plateau starts (see Fig. 2) and thus where the 
merging events can take place. As the shock merging dis- 
turbance occurs in a progressively thinner atmosphere it 
generates weaker and weaker resonance oscillations, which 
in turn decreases the number of shock merging events with 
time. 

Figure 3 depicts the shock merging heights as function 
of time for two adiabatic wave computations with P = 
15 s (diamonds) and P = 30 s (plus signs). As in the 30 s 
case the wavelength is much larger than in the 15 s case, 
it is seen that only 4 shock merging events happen in our 
computational domain, which extends to 2220 km, and 
even then, shock merging occurs only in the initial switch- 
on phase at times t < 700 s. This allows the resonances 
to die out completely. In the 15 s case there are 108 shock 
merging events, but it is seen that with time these merging 
events occur at increasingly greater height. As discussed 
before, there will be a time when shock merging no longer 
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Fig. 4. Power spectra (left and middle columns) and velocities at height z = 2000 km (right column) for an excitation with 
monochromatic waves of period P = 15, 20, 30 s (from top to bottom) and flux F\ = 1 ■ 10 8 erg cm -2 s" 1 , for non-adiabatic 
waves. The power spectra at height z = km, scaled by 0.1, are shown dotted, those at height z = 2000 km drawn. The velocity 
is sampled from 1500 to 3548 s in 1 s intervals. 



occurs in the computational domain. Thus in the 15 s case 
the resonances will die out as well. 

In summary we conclude that for adiabatic monochro- 
matic waves, smaller wave periods lead to more shock 
merging events which generate resonance oscillations. But 
the unrestricted heating in the adiabatic case modifies the 
atmospheric structure in such a way that eventually, at 
a given height, further shock merging is prevented and 
the resonance oscillations die out. Thus for adiabatic non- 
linear monochromatic wave excitations, regardless of the 



wave period, the resonances will eventually die out and 
the forced oscillations will prevail. 

3.3. Monochromatic non-adiabatic wave excitation 

We now consider wave calculations with radiation damp- 
ing. Figure 4 (left two columns) depicts the power spectra 
at z = 0, 1300 and 2000 km for waves of period P = 15, 20 
and 30 s. These spectra differ in two important ways from 
those of Fig. 2 of Paper II. They are now computed allow- 
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ing for radiation damping and the time span of the Fourier 
analysis is now from 1500 to 3548 s instead of from 500 to 
2048 s. The velocity amplitudes of the waves, displayed in 
the right column of Fig. 4 as well as in the corresponding 
Fig. 2 of Paper II, show that resonance oscillations occur 
well past the decay time of the initial switch-on effect. 

Comparison of the two P — 15 s spectra at z — 
2000 km shows that the adiabatic calculation in Fig. 1 
has a much sharper resonance peak than that of Fig. 4. 
This is because in Fig. 1, the Fourier analysis is over the 
time span 500 to 2548 s, which still includes a considerable 
portion of the switch-on effect, while in Fig. 4 in the time 
between 1500 to 3548 s this effect is negligible. Another 
reason is that radiation damping, by decreasing the wave 
amplitudes, decreases the shock speeds and thus counter- 
acts shock merging. Less shock merging always implies a 
larger amount of high-frequency power. 

3.4- Effect of the generated chromosphere model, non- 
adiabatic case 

Including radiation changes the behaviour of the chro- 
mosphere model drastically. The temperature increase in- 
duced by shock heating is now counteracted by radiative 
cooling such that after some time, heating and cooling bal- 
ance each other. Figure 5 shows the mean temperatures as 
function of height at 1500 s intervals both for the P = 15 
and 30 s waves. It is seen that in both cases after about 
2000 s, a dynamically stable mean chromosphcric temper- 
ature rise is established up to a height of z = 1400 km. 
Here shock heating is balanced by H~ and Mg II k line 
cooling. At greater heights, where Mg II cooling becomes 
inefficient and Lya cooling becomes appreciable, a sta- 
ble dynamical temperature distribution is achieved after 
about t = 5000 s. 

It should be noted, however, that because of various 
reasons we do not consider this established temperature 
profile in the high chromosphere to be fully realistic. First 
we have used an arbitrarily selected monochromatic wave, 
second our treatment of radiation losses is rather incom- 
plete both because of technical reasons and the small num- 
ber of coolants considered. In addition, we have neglected 
the ionization energy of hydrogen in the energy budget 
of the wave and have omitted a fully consistent time- 
dependent treatment of the H ionization and recombina- 
tion. In our treatment the ionization equilibrium follows 
the temperature variations instantaneously due to our so- 
lution of the time-dependent statistical rate equations for 
hydrogen. As shown by Carlsson & Stein (1992, 1994), 
the incorporation of the above mentioned effects signifi- 
cantly influences the temperature profiles particularly in 
the upper solar chromosphere. 

Let us now consider the acoustic wave spectra in Fig. 4 
in view of the dynamic chromospheric models obtained. 
Figure 4 shows at height z = 2000 km that the resonance 
oscillations in all three wave cases still persist after t = 
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Fig. 6. Shock merging heights as function of time for radia- 
tively damped monochromatic waves with an acoustic energy 
flux of Fa = 1 ■ 10 8 erg cm~ 2 s _1 and a period of P — 15 s 
(diamonds) as well as P — 30 s (plus signs) . 



3500 s, whereas at the height z = 1300 km, the spectra 
consist primarily of monochromatic components. Figure 6 
gives the height of the shock merging events as function 
of time for the P = 15 and 30 s waves. Comparison with 
Fig. 3 shows that, while the rate of shock merging for the 
30 s wave decreases with time similarly as in Fig. 3, that of 
the 15 s wave is very different. In the latter case the rate of 
shock merging events above z = 1400 km does not appear 
to decrease with time. This shows that in the cases P = 15 
and 20 s, a strong resonance contribution with very little 
or no contribution of the initial monochromatic excitation 
becomes established which persists with time. 

As the shock merging events for the P = 30 s wave oc- 
cur increasingly close to the upper computational bound- 
ary at z = 2220 km it can be expected that for this 
wave the spectrum at z = 2000 km in Fig. 4 will lose 
its resonance contributions and eventually show only 
the monochromatic contribution. This supports the con- 
cept already discovered in Paper II of a critical period 
P cr « 25 s below which the resonance oscillations are self- 
sustaining, while above that period the resonances die out. 

To summarize we find that for non-adiabatic (i.e., ra- 
diatively damped) monochromatic waves the resonance 
behaviour is different compared to that of the adiabatic 
monochromatic waves. As the radiative cooling allows 
the atmosphere to reach a dynamical steady state, shock 
merging processes take place now in an atmosphere of a 
fixed mean temperature. Here waves with periods P < P cr 
find the sound speeds low enough to catch up with one an- 
other in an atmospheric oscillation kicked on by the initial 
switch-on effect. These shock merging events are able to 
perpetuate the resonance oscillations indefinitely. For pe- 
riods P > P cr , shock merging becomes increasingly rare 
allowing the resonances to die out. 
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Fig. 5. Mean temperatures as function of height at times indicated 1 = 500 s, 2 = 2000 s, 3 = 3500 s, 4 = 5000 s for radiatively 
damped monochromatic waves with a wave energy flux of Fa = 1 ■ 10 8 erg cm~ 2 s" 1 and a period of P = 15 s (left panel) as 
well as P — 30 s (right panel). The temperature of the initial atmosphere is shown dotted. 



3.5. Excitation by acoustic spectra for non-adiabatic waves 

We now consider the excitation of the solar atmosphere 
by various acoustic wave spectra with radiation damping 
included. For this exploratory investigation we select only 
two of the four excitation spectra considered in Paper II. 
The first one is a Gaussian spectrum with a central peak at 
frequency uic = 27r/35 Hz, extending from 10 to 50 mHz, 
the second is a stochastic spectrum extending from 5 to 
55 mHz (see Paper II) . As in Paper II we represent the 
input spectrum by 101 partial waves equidistantly spaced 
in frequency with Av = 0.5 mHz. The lowest frequency 
point is about a factor of 0.7 below the cut-off frequency 
va ~ 5 mHz, the highest about a factor of 11 larger. 

For a comparison with our monochromatic studies we 
take the same total energy flux Fa — 1 • 10 8 erg cm~ 2 s _1 
for the spectra. The initial acoustic wave spectra which 
excite the solar atmosphere model at height z = km are 
shown (dotted) in all panels of Figs. 7 and 8. These fig- 
ures show (drawn) the acoustic wave spectra at heights 
z = 1000, 1300 and 2000 km (from top to bottom). Fur- 
thermore, different timespans for the Fourier analysis are 
taken. We selected the time intervals 500 - 2548 s (left 
column) and 2500 - 4548 s (right column). There are two 
striking results from these spectra. First, it is seen that in 
both wave cases there is a strong and with height progress- 
ing shift of the spectral power from the initially specified 
one to one which shows a predominant contribution in 
the 6-8 mHz (3 min) range. Comparing Figs. 7 and 8 it 
is seen that despite the very different initial shapes of the 
two spectral cases, they become increasingly similar for 
greater height. At 2000 km, one has almost pure 3 min 
spectra in both cases. 

The second striking feature of Figs. 7 and 8 is that 
the spectra at a given height apparently do not vary with 
time. Comparing the spectra computed from the time in- 



terval 500 - 2548 s with those of 2500 - 4548 s at a given 
height shows that there is essentially no change. This is 
very different from the monochromatic adiabatic results 
where the resonance oscillations are transient and die out 
with time. 

As found in our monochromatic studies above, a com- 
bination of two conditions will lead to a continuous pres- 
ence of resonance oscillations: the atmosphere must reach 
a dynamical equilibrium and the shock merging events 
must go on indefinitely. Both conditions are met in an ex- 
citation with frequency spectra. Figure 9 shows the mean 
temperatures of the atmosphere as function of time for the 
two wave cases considered. It is seen that like in Fig. 5, 
a dynamical steady state temperature distribution is es- 
tablished below 1400 km height after 2000 s and some- 
what later after 3500 s at greater height. For the second 
condition, Fig. 10 shows the plots of the shock merging 
heights of the two wave cases considered. It is seen that for 
the Gaussian spectrum, shock merging goes on indefinitely 
with time in the height interval between 1000 km and the 
computational boundary at 2220 km. For the stochastic 
wave case, incessant shock merging occurs in the range 
600 to 2000 km. 

Because shock merging implies the presence of shocks, 
and because the initial acoustic frequency spectra have 
relatively small amplitude with no shocks present, it is 
clear that the shocks have to be generated during the wave 
propagation. Figure 11 shows for the stochastic case that 
different to the situation of monochromatic waves, where 
all shocks form close to 400 km height (see Fig. 2), the 
shocks in the case of an acoustic spectrum are formed in a 
wide range between 300 and 1100 km height. Above 1100 
km height, the waves are all sawtooth-type shock waves. 
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Fig. 7. Power spectra at (dotted), 1000, 1300, 2000 km height (from top to bottom), Fourier analysing the time span 500 - 
2548 s (left column) and 2500 - 4548 s (right column) for an atmosphere excited at height z = km by a Gaussian acoustic 
frequency spectrum with a flux of Fa = 1 • 10 s erg cm -2 s _1 and a central maximum at period P = 35 s. The wave calculation 
is non-adiabatic. 



3.6. Shock merging and the generation of 3 min oscilla- 
tions 

At this point we want to discuss in more detail the rela- 
tion between shock merging and the generation of 3 min 
oscillations. In the work of Lites et al. (1993) as well as 
Cheng & Yi (1996) it is argued on basis of observations 
that in the solar atmosphere the 3 min signal is already 
present in the lower photosphere and that therefore shock 
merging is not an important process for the generation of 
3 min oscillations. We believe that contrary to this view, 



our present calculations clearly show that shock merging 
is a powerful method for generating 3 min oscillations. 

Our monochromatic calculations indicate that if the 
wave periods are not too small, the resonance oscilla- 
tions die out. However, our calculations with acoustic 
wave spectra show a fundamentally different behaviour. 
One might think that in an acoustic spectrum (cf. Fig. 2) 
the individual partial waves simply form a wavetrain of 
sawtooth waves and thus do not generate low frequency 
oscillations. However, the stochastic nature of the super- 
position of partial waves in an acoustic spectrum perpetu- 
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Fig. 8. Same as Fig. 7, however for a stochastic acoustic frequency spectrum. 



ally leads to large amplitude fluctuations which invariably 
generate shocks of different wave period and strength. In 
addition, note that stochastic superpositions lead to per- 
petual resonances even in the small amplitude case (see 
Fig. 14 of Paper I). 

The propagation of the shocks in an oscillating atmo- 
sphere subsequently results in shock merging events which 
in turn drive the atmospheric oscillations. Thus acoustic 
wave spectra will invariably lead to perpetual resonance 
oscillations. As in our assumed input spectra the 3 min sig- 
nal is absent, these resonance oscillations must have been 
generated by the atmosphere in response to the wave exci- 
tation. As the acoustic spectrum calculations of Musielak 
et al. (1994) do not show a prominent 3 min contribution 



at the top of the convection zone, we consequently con- 
clude that the 3 min band present in the Lites et al. ob- 
servations must have been added later by hydrodynamic 
effects in the solar atmosphere. This conjecture will be 
investigated further in the subsequent paper V. 

3. 7. Acoustic spectra and the dynamic structure of the at- 
mosphere 

A very important result found by comparing the two cases 
of Fig. 9 is that the mean temperature of the atmosphere 
in a dynamical steady state strongly depends on the initial 
wave spectrum. This result also applies to the temperature 
behaviour of Fig. 5 for the unrealistic case of monochro- 
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Fig. 10. Shock merging heights as function of time for the Gaussian acoustic wave spectrum (left panel) and the stochastic 
spectrum (right panel). 



matic waves. We therefore conclude that the structure of 
stellar chromospheres is intimately determined by the ini- 
tial acoustic wave spectrum, generated in the convection 
zone. It thus should be possible to infer the initial wave 
spectrum by detailed simulations by comparing theoreti- 
cal spectra (sec e.g. Figs. 7 and 8) and observed spectra 
at various atmospheric heights. 

4. Conclusions 

From our adiabatic and non-adiabatic wave calculations 
for a solar atmosphere model, excited at the bottom by 
large amplitude monochromatic waves as well as by acous- 
tic frequency spectra, employing a realistic mechanical 
flux of Fa = 1 • 10 8 erg cm -2 s _1 , we draw the following 
conclusions: 

1. Adiabatic wave calculations, because of the unbalanced 
heating, invariably lead to chromospheric temperature 



plateaus with perpetually rising mean temperatures. This 
time-dependent growth of the mean temperature makes 
it increasingly difficult for shock merging to occur. The 
consequence is that the 3 min type resonances, kicked on 
by the shock merging events, die out and after some time 
only the monochromatic signal survives. 

2. Non-adiabatic calculations, where the radiative cooling 
by NLTE H~ continua as well as by Mg II k and H 
Lya lines are considered, invariably establish a dynam- 
ically generated mean chromospheric temperature distri- 
bution, which at heights below 1400 km is established after 
2000 s and at greater height after about 3500 s. 

3. For excitation by non-adiabatic monochromatic waves, 
as already noted in Paper II, a critical frequency v cr ~ 
1/25 Hz is found, which separates domains of drastically 
different resonance behaviour. For frequencies v < u cr , 
the atmospheric resonance decays, leaving behind only the 
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Fig. 11. Shock formation heights as function of time for a 
radiatively damped stochastic wave spectrum of an acoustic 
energy flux of Fa = 1 • 10 8 erg cm~ 2 s _1 . 



monochromatic signal, while for v > v cr the resonance os- 
cillation persists and is kicked on by shock merging. For 
wave periods P < 25 s, the shock merging occurs undi- 
minished in the height range 1400 to 2200 km. Here the 
resonances appear self-sustaining. 

4. Using two different acoustic wave spectra (Gaussian 
and stochastic) to excite the atmosphere at height z = 
km, we found that in non-adiabatic wave calculations, 
the spectra showed a characteristic shift towards lower 
frequency such that at height z = 2000 km essentially 
only a pure 3 min band exists. 

5. We also found that at any given height the acoustic 
spectra, after an initial time of 500 s, did no longer change 
with time. 

6. In the height interval 1000 - 2200 km for the Gaus- 
sian and 600 - 2000 km for the stochastic spectrum, shock 
merging occurs persistently with time. This self-sustaining 
shock-merging behaviour is responsible for the fact that 
there is a 3 min resonance band in the acoustic spectra, 
which becomes increasingly pronounced with atmospheric 
height. 

7. From our calculations we surmise that the 3 min com- 
ponent observed by Lites et al. (1993) in a low lying iron 
line, which is not found in the acoustic spectrum gener- 
ated in the convection zone (Musielak et al. 1994), has 
been added later by the solar atmosphere as response to 
the propagation of the acoustic wave spectrum. 

8. We find that the dynamically generated mean chromo- 
sphcric temperature structure is strongly dependent on 
the assumed initial acoustic wave spectrum. This indicates 
that detailed wave simulations and height-dependent 
spectral observations will allow to empirically determine 
the velocity fluctuations in the acoustic wave generation 
region. This could be an important independent check of 
existing convection and sound generation simulations. 



9. The appreciable atmospheric expansions found in Paper 
II did not occur in the present simulations because we now 
use an Eulerian code. 
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